NONLINEAR DYNAMICS OF HIGH-BRIGHTNESS BEAMS 
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Abstract 



The consideration of transverse dynamics of relativis- 
tic space-charge dominated beams and halo growth due to 
bunch oscillations is based on variational approach to ratio- 
nal (in dynamical variables) approximation for rms enve- 
lope equations. It allows to control contribution from each 
scale of underlying multiscales and represent solutions via 
exact nonlinear eigenmodes expansions. Our approach is 
based on methods provided possibility to work with well- 
localized bases in phase space and good convergence prop- 
erties of the corresponding expansions. 

1 INTRODUCTION 

In this paper we consider the applications of a new nu- 
merical-analytical technique based on the methods of local 
nonlinear harmonic analysis or wavelet analysis to nonlin- 
ear rms envelope dynamics problems which can be char- 
acterized by collective type behaviour [1], [2]. Such ap- 
proach may be useful in all models in which it is possible 
and reasonable to reduce all complicated problems related 
with statistical distributions to the problems described by 
systems of nonlinear ordinary/partial differential equations 
with or without some (functional) constraints. Wavelet 
analysis is a set of mathematical methods, which gives the 
possibility to work with well-localized bases in functional 
spaces and gives the maximum sparse forms for the general 
type of operators (differential, integral, pseudodifferential) 
in such bases. Our approach is based on the variational- 
wavelet approach from [3]-[14], that allows to consider ra- 
tional type of nonlinearities in rms dynamical equations. 
The solution has the multiscale/multiresolution decompo- 
sition via nonlinear high-localized eigenmodes, which cor- 
responds to the full multiresolution expansion in all under- 
lying time/space scales. We may move from coarse scales 
of resolution to the finest one for obtaining more detailed 
information about our dynamical process. In this way we 
give contribution to our full solution from each scale of 
resolution or each time/space scale or from each nonlin- 
ear eigenmode. The same is correct for the contribution 
to power spectral density (energy spectrum): we can take 
into account contributions from each level/scale of reso- 
lution. Starting in part 2 from general rms envelope dy- 
namics model we consider in part 3 the approach based on 
variational-wavelet formulation. We give explicit represen- 
tation for all dynamical variables in the base of compactly 
supported wavelets or nonlinear eigenmodes. Our solutions 
are parametrized by the solutions of a number of reduced 
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algebraical problems one from which is nonlinear with the 
same degree of nonlinearity and the others are the linear 
problems depend on particular wavelet type approach. In 
part 4 we consider numerical modelling based on our ana- 
lytical approach. 

2 RMS EQUATIONS 

We consider an approach based on the second moments 
of the distribution functions for the calculation of evolution 
of rms envelope of a beam. The rms envelope equations are 
the most useful for analysis of the beam self-forces (space- 
charge) effects and also allow to consider both transverse 
and longitudinal dynamics of space-charge-dominated rel- 
ativistic high-brightness axisymmetric/asymmetric beams, 
which under short laser pulse-driven radio-frequency pho- 
toinjectors have fast transition from nonrelativistic to rela- 
tivistic regime [1]. The analysis of halo growth in beams, 
appeared as result of bunch oscillations in the particle-core 
model, is also based on three-dimensional envelope equa- 
tions [1], [2]. Let f{xi) be the distribution function, which 
gives full information about noninteracting ensemble of 
beam particles regarding to trace space or transverse phase 
coordinates {xi, Xj). Then we may extract the first nontriv- 
ial effects of collective dynamics from the second moments 



RMS emittances are given by 
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We consider the following most general case of rms enve- 
lope equations, which describe evolution of the moments 
in the model of halo formation by bunch oscillations (ref. 
[2] for full designation): 
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where (Jx{s),ay{s),az{s) are bunch envelopes, ^£!;,<^j,, 
^z = F{ax,a-y,o-z)- After transformations to Cauchy form 
we can see that all these equations from the formal point of 
view are not more than ordinary differential equations with 
rational nonlinearities and variable coefficients. Also, we 
consider regimes in which we are interested in constraints 
on emittances: 

(4) 
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where a are constants. In the same way according to 
[2] we may consider the case of energy-type functional- 
differential constraints on emittances. A different approach 
is considered in our related paper in this Proceedings [15]. 

3 RATIONAL DYNAMICS WITH 
CONSTRAINTS 

Our problems above may be formulated as the systems 
of ordinary differential equations 
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with initial (or boundary) conditions Xi{Q), Xi{T) and 
Pi , Qi are not more than polynomial functions of dy- 
namical variables Xj and have arbitrary dependence on 
time/length parameter. Of course, we consider such Qi (x) 
which do not lead to the singular problem with Qi{x), 
when t ^ OT t ^ T, i.e. Q,{x{0)), Q,{x{T)) ^ oo, 
0. We'll consider these equations as the following oper- 
ator equation. Let L be an arbitrary nonUnear (rational) 
matrix differential operator of the first order with matrix 
dimension d (d=6 in our case) corresponding to the sys- 
tem of equations (5), which acts on some set of functions 
* = = t e n C i? from 

= = 0, (6) 

where R = R{t, d/dt, ^f). 

Let us consider now the N mode approximation for so- 
lution as the following ansatz (in the same way we may 
consider different ansatzes); 
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We shall determine the coefficients of expansion from the 
following variational conditions (different related varia- 
tional approaches are considered in [3]-[14]): 
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We have exactly dN algebraical equations for dN un- 
knowns Ur- So, variational approach reduced the initial 
problem (5) or (6) to the problem of solution of functional 
equations at the first stage and some algebraical problems 
at the second stage. Here ip^ (t) are useful basis functions of 
some functional space {L^,LP, Sobolev, etc) correspond- 
ing to concrete problem. As result we have the following 
reduced algebraical system of equations (RSAE) on the set 
of unknown coefficients of expansions (7): 
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(9) 



possible time dependence) of LHS of initial system of dif- 
ferential equations (^ and as consequence are coefficients 
of RSAE. Pij are the coefficients (with possible time de- 
pendence) of RHS of initial system of differential equa- 
tions (^ and as consequence are the coefficients of RSAE. 
/ = (ii, ...,iq+2), J = (ji, jp+i) are multiindexes, 
by which are labelled ai and f3j, the other coefficients of 
RSAE (|): 

i<ifc<p+i 

where p is the degree of polynomial operator P (|]) 
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where q is the degree of polynomial operator Q (|^), ie = 

(l,...,(7 + 2), ^di^^Jdt. 

According to [3]-[14] we may extend our approach to 
the case when we have additional constraints (4) on the set 
of our dynamical variables or x. In this case by using 
the method of Lagrangian multipliers we again may apply 
the same approach but for the extended set of variables. As 
result we receive the expanded system of algebraical equa- 
tions analogous to the system (9). Then, after reduction we 
again can extract from its solution the coefficients of ex- 
pansion (7). Now, when we solve RSAE (^) and determine 
unknown coefficients from formal expansion (7) we there- 
fore obtain the solution of our initial problem. It should 
be noted if we consider only truncated expansion (7) with 
N terms then we have from (^ the system of x d al- 
gebraical equations with degree £ = max{p, q} and the 
degree of this algebraical system coincides with degree of 
initial differential system. So, we have the solution of the 
initial nonlinear (rational) problem in the form 
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x{t)=x{0)+J2''kMt), 
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where coefficients are roots of the corresponding re- 
duced algebraical (polynomial) problem RSAE (^. Con- 
sequently, we have a parametrization of solution of initial 
problem by solution of reduced algebraical proble m (^ . 

The problem of computations of coefficients aj (|ll|), /3j 
diol ) of reduced algebraical system may be explicitly solved 
in wavelet approach. The obtained solutions are given in 
the form (|l2), where ipkit) are wavelet basis functions. In 
our case ipk (t) are obtained via multiresolution expansions 
and represented by compactly supported wavelets. Because 
affine group of translation and dilations is inside the ap- 
proach, this method resembles the action of a microscope. 
We have contribution to final result from each scale of res- 
olution from the whole infinite scale of spaces: 

...V-2 C V^i C C C ^2 C 



where operators L and M are algebraization of RHS and 
LHS of initial problem (^). Qij are the coefficients (with 



where the closed subspace Vj{j G Z) corresponds to level 
j of resolution, or to scale j. This multiresolution functional 



space decomposition corresponds to exact nonlinear eigen- 
mode decompositions (12). 

It should be noted that such representations give the 
best possible localization properties in the corresponding 
(phase)space/time coordinates. In contrast with differ- 
ent approaches formulae (7), (12) do not use perturba- 
tion technique or linearization procedures and represent 
dynamics via generalized nonlinear localized eigenmodes 
expansion. So, by using wavelet bases with their good 
(phase)space/time localization properties we can construct 
high-localized (coherent) structures in spatially-extended 
stochastic systems with collective behaviour. 
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Figure 1: Multiscale decomposition. 




Figure 2: ax-<Jy section. 



4 MODELLING 



So, our N mode construction (7), (12) gives the follow- 
ing representation for solution of rms equations (3): 



xt^{t) + J2^'i^'^t)^ (13) 

i>N 



where x^{s) may be represented by some family of (non- 
linear) eigenmodes and gives as a result the multiresolu- 
tion/multiscale representation in the high-localized wavelet 
bases. The corresponding decomposition is presented on 
Fig. 1 and two-dimensional transverse section ax-dy on 
Fig. 2. 
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